Time of avalanche mixing of granular materials in a half filled rotated drum 
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The avalanche mixing of granular solids in a slowly rotated 2D upright drum is studied. We 
demonstrate that the account of the difference S between the angle of marginal stability and the 
angle of repose of the granular material leads to a restricted value of the mixing time r for a half 
filled drum. The process of mixing is described by a linear discrete difference equation. We show 
that the mixing looks like linear diffusion of fractions with the diffusion coefficient vanishing when 
5 is an integer part of n. Introduction of fluctuations of 5 supresses the singularities of t{5) and 
smoothes the dependence t{S). 
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The evolution of granular materials in rotated drums is 
one of the most interesting topics of soft matter physics 
Granular flows induced in such a way demon- 
strate a lot of intriguing phenomena including, e.g., seg- 
regation of different fractions, mixing, self-organization, 
self-organized criticality, etc [p"2|-^. In principle, the 
system is three-dimensional, and the process features 
depend strongly on the particular characteristics of the 
grains involved. A theoretical description of these pro- 
cesses is still unsolved problem, and only the simplest 
cases were considered analytically. 

Decrease of the dimensionality of the system by one 
simplifies the consideration, although docs not make the 
arising problems to be less interesting. The avalanche 
mixing phenomenon fl^ observed in a slowly rotated 2D 
upright drum, i.e., in a flat disc, turned to be so impres- 
sive that one of figures of the paper has made a front 
page of the Nature. 

Let us introduce the avalanche mixing problem 

m — 



^Jl-25|. A partially filled 2D upright drum rotates 
very slowly around its axis (see Fig. 0). Small grains of 
the granular material move as a whole with the drum. 
There is, however, one exception: periodically, when the 
grains appear on the free surface, they fall down along 
it. If, initially, there is several separated fractions of the 
material, the grains of the different fractions may mix 
only in the avalanches which fall down along the free 
surface. One has to describe the evolution of the mix- 
ing, i.e., to describe how the process develops. The main 
quantity under study is the time of mixing r, that is the 
characteristic time after that the material will be mixed 
homogeneously. It may be introduced, since the relax- 
ation to the final homogeneous state in such a process 
proceeds by an exponentional law. 

There are two particular values of the angle of the 
free surface slope of a granular material - the angle of 
marginal stability and the angle of repose. (We assume 
that the free surface is straight. That is reasonable for 
a slowly rotated drum.) When the free surface slope 



achieves the angle of marginal stability, the grains from 
the sector of the angle S near the upper half of a free 
surface fall down (see Fig. |l| - grains fall down from sec- 
tor u to sector d). Here, S is the difference between the 
angle of marginal stability and the angle of repose. The 
avalanche stops when the free surface slope approaches 
the repose angle, and the drum has to turn by the angle 
S to launche the next avalanche. Therefore, the only rea- 
son for mixing in the problem under consideration is the 
granular flow (avalanches) between these wedges. 

We assume that the fractions become to be mixed ho- 
mogeneously after the fall of the avalanche - the mixture 
in sector d in Fig. |l|. This case, in which the avalanche 
mixing proceeds most quickly, appears to be near the ex- 
perimentaly realized situation fl^ . At least, using this 
assumption, one may obtain the minimal possible time 
of the avalanche mixing. 

When (5 — > 0, the mixing was described in the frames 
of a simple geometrical approach [ p2pl|j2^ , p^ , p5[ , and an 
analytical theory |2^j2^ describing the experimental ob- 
servations [|l2| surprisingly well has been proposed. Nev- 
ertheless, there is a very interesting particular case of the 
half filling, in which one certainly can not assume that S 
is zero. (Note, that in the experiment of Metcalfe et al 
|l2| , (5 « 8°.) Indeed, this assumption leads immediately 
to an infinite value of the mixing time for a half filled 
drum: there is no any mixing between the grains of the 
adjacent wedges ii S — 0, and the material will never mix. 

The finite difference between the angle of marginal sta- 
bility and the repose angle makes r to be restricted even 
if the drum is half filled. Hence, the finiteness of 6 de- 
termines the avalanche mixing process in this particular 
case that is considered here. 

Let us assume that there are two fractions - black and 
white, and that the small grains of the different fractions 
are distinguished from each other only by their color. 
The material in sector d of the angle S is mixed homo- 
geneously after the avalanches falls. Then, after the first 
half of turn, the mixing dynamics can be described by a 
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set of concentration values, {cm}, where Cm is the con- 
centration of the black fraction in sector d (Fig. ^ after 
the m-th avalanche falls, that happens at time t — mS. 
Here, time is an angle of drum rotation. Note that a 
drum radius is not essential for the consideration and 
does not appear in our expressions. In the case of a half 
filled drum, the equation for c„i turns to be very simple: 

(1) 

It can be understood easily from Fig. |l|. The grains 
of only two adjacent wedges mix. In each of these 
wedges the material is homogeneous after the previous 
avalanches, so two terms in the right part of Eq. (|^) 
are the contributions of the black grains from one of 
these two wedges, taken with the corresponding weight 
~ their relative angles in sector u in Fig. |^. Here, 
[ ] denotes an integer part of a number. The values 
Cm , m — — [tt/J] , — [tt/J] + 1 1 , may be taken as 
an initial condition if one does not consider the first half 
of turn. Obviously, the initial conditions are not essential 
for the mixing time value. 

For brevity, we use the notation n = [tt /S\ for the in- 
teger part and a = tt / 5 — [t^/S\ for the fractional part. 
With these notations, Eq. (]l|) may be written in a more 
compact form: 



(2) 



According to usual prescriptions j2^], one may search 
for the solution of Eq. (|l|) or Eq. in the form of a lin- 
ear combination of the terms A™, where \j are the roots 
of the characteristic polynomial 



(1 -a)A 







(3) 



To obtain the long time relaxation we should find, among 
the roots Aj of Eq. (^, the maximal value |Aj| = Xmax 
that is lower then the root A = 1 (note that there is no 
multiple roots; A = 1 corresponds to a stationary situ- 
ation). The mixing time r is expressed in the terms of 
Amaa;: = (1/^ ) log(l/ Amaa; ) , SO it may bc Calculated 
directly for any value of a and a finite n. The result is 
shown by the solid line in Fig. ||. One can see that the 
mixing time is infinite for only integer values of tt/J. It 
has local minima near the points tt/J = k -\- 1/2, where 
/J = 2,3,.... 

The process of the mixing may be interpreted in the 
following way. It follows from Eq. (||) that the long scale 
and the short one are separated if n ^ 1. Therefore, in 
such a situation, it is reasonable to use different variables 
for these two scales. Let us introduce a new time vari- 
able t. We observe periodically the system under study 
with the interval of periodicity tt (e.g., using flashes like 
in the stroboscopic effect) and describe the result of the 



observation by the quantity c((/?,f), the concentration of 
the black fraction at the angle ip at the "time" t. Here, 

< (/? < TT. 

At n ^ 1, in the continuous limit, c„j — > c{m5). Sub- 
tracting c((m — a)(5) from both sides of Eq. (H) and using 
the obvious relation, {n + a)5 — tt, we obtain immediately 
the equation for c((^,t): 



9c((^,t) 5^a{l-a) d'^c{ip,i) 



dt 
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(4) 



Hence, in such an approach, the mixing looks like the 
diffusion governed by Eq. (^). We do not write out its 
solution since only the mixing time, i.e., the time of the 
exponentional relaxation to a homogeneous state, is in- 
teresting for us. It follows from the diffusion coefficient 
and is of the form: 



(5) 



This result (which is obtained in the limit (5 <C tt) is 
shown by the dashed curve in Fig. 2. One may see that 
the deviations from the results obtained directly from 
the characteristic polinomial equation are small even for 
large 5. The considerable difference is visible only in the 
unphysical range, 7r/3 < 5 < n/2. 

Above, we considered the fixed S. Nevertheless, really, 
in an experiment, the angle S is is not fixed but fiuctuates 
with time. A particular form of the distribution function 
of i5 is dependent on characteristics of the mixed grains 
and is not discussed here. What is the effect of these 
ffuct nations? 

Let the difference 6m between the marginal stability 
angle and the repose angle for the m-th avalanche be 
fiuctuating and, hence, depending on m. Let the dis- 
tribution function for it be P{S). One introduces the 
total angle of the drum turn, Om — '^iLi ^i- Then, after 
the m-th avalanche, the wedge between 9m-i and 9m of 
width 5m is homogeneously mixed with the concentration 
of the black fraction Cm- 

The material for the mixing is taken from the wedge 
Om-i — TT < t < 9m — T^- Now, uulike the previ- 
ously considered case without the fluctuations, it may 
be more than two wedges with different concentration of 
the black grains in this sector. Therefore, if we assume 



7 m- 1 



< 9j+k < 



TT, where the value 



of k depends on the particular set of {<5m}, one obtaines 
the following equation: 

Cm = ^ {(^rn-l - 7r)]Cj + Sj+iCj + i + . . . + 

Sj+kCj+k + [{9m - tt) - 9j+k]cj+k+i} (6) 

that is the strict generalization of Eq. (^. We study this 
equation numerically with account of the fluctuations of 
dm to find the mixing time. The results for the homoge- 
neous distributions P{6) = e{S - {1 - A)S)e{{l + A)S - 
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(5)/ (2 A) with A = 1/8 and A = 1/4 are presented in 
Fig. H (here, 6( ) is the theta- function) . (Note, that, 
in Fig. H, we do not connect the points obtained from 
the numerics.) One sees that the singularities of r are 
effectively suppressed by the fluctuations. The oscilla- 
tions of the dependence decrease with decrease of 6, and 
{tt / {26'^))t~^ approaches a constant value at small 6. For 
small values of the parameter A, this value is approached 
more quickly than for large ones but it is easy to check 
that, for all C > 0, it equals approximately 0.175. That 
is close to 1/6 — daa{l — a), i.e., to the average value 
of (7r/(2(52))T-i at smah S. 

In the experiment |jl^, S — 8° ± 2°. From our esti- 
mation, that leads to r ~ 70 turns of the drum. There 
are no available experimental data yet to be compared 
with this value. The measuring of r in such a situation 
would be really interesting. One can expect that our esti- 
mation is quite reasonable because the agreement of the 
results of the geometrical approach for (5 = with 
the experiment is excellent. 

In conclusion, we have described analytically the dy- 
namics of the avalanche mixing in the case of a half 
filled slowly rotated drum when the difference d between 
the angle of the marginal stability and the angle of re- 
pose plays the principal role determining the value of the 
mixing time. The mixing looks like diffusion of grains 
between wedges with the diffusion coefficient vanishing 
when the angle S is an integer part of tt. In these points, 
the dependence r((5) has singularities which may be effec- 
tively suppressed by fluctuations of S. The proposed ap- 
proach is, in fact, geometrical, and does not contain any 
other parameters apart of S or its distribution. Two main 
assumptions were made: (i) We considered the granular 
material consisting of very small grains (much smaller 
than the drum size), (ii) The mixing of the fractions 
after each elementary avalanche were assumed to be ho- 
mogeneous. These assumptions reserve space for future 
study of the avalanche mixing problem. 
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FIG. 1. Schematic illustration of the avalanche mixing in 



a half filled 2D drum. Mixing begins from an inhomogeneous 
distribution of the fractions, e.g., from the configuration (A). 
(B) The drum starts rotate anticlockwise. When the free sur- 
face slope attains the marginal stability angle value, the grains 
of different fractions flow quickly from sector u to sector d, 
undergoing mixing. Then the free surface slope turns to be 
equal to the repose angle of the material. The angle 5 is the 
difference between the angle of marginal stability and the an- 
gle of repose. In the figure, the angle of repose is zero. If 5 is 
not an integer part of tt or it fluctuates with time, finally the 
material is homogeneously mixed (C). 
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FIG. 2. Scaled inverse avalanche mixing time {n/25 )t~^ 
Ysn/5. The solid line describing the case without fluctuations 
of S is obtained from the solution of the characteristic poly- 
nomial, Eq. (^). The dashed line represents Eq. (P) which 
is, strictly speaking, found for small 5 <^ n. The points are 
results of the direct simulation of the evolution of the mix- 
ing with account for the fluctuations of S, Eq. was used. 
The squares correspond to the homogeneous distribution of S 
between (1 — 1/8)5 and (1-1-1/8)5, the circles - to the homo- 
geneous distribution of S between (1 — 1/4)5 and (1+1/4)5. 
The lines between the points are not shown. 
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